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ON  THE  RESPONSE  OF  KNOLLENBERG  AEROSOL  COUNTERS 

by 

R.  R.  Allan 
P.  C.  Ashdown 


SUMMARY 

It  is  known  that  the  response  of  light-scattering  aerosol  counters  is 
sensitive  to  the  aerosol  refractive  index,  and  that  even  for  a  given  refrac¬ 
tive  index  there  are  ’multi-valued’  regions  where  significantly  different 
sizes  of  particle  give  the  same  response.  In  natural  aerosols,  the  particle 
size  distribution  inferred  directly  from  measurements  shows  a  'knee' 
followed  by  a  sharp  drop.  Moreover  the  p  >sition  of  the  knee  should  depend 
on  the  aerosol  refractive  index. 


In  practice  the  most  obvious  feature  is  a  knee  near  1  um  nominal  size 
on  range  O  of  the  Knollenberg  ASAS  counter,  although  further  features  are 
present  at  much  poorer  resolution  on  range  1  of  the  Knollenberg  CSAS 
counter.  Here  we  show  that  the  position  of  the  knee  correlates  very  well 
with  relative  humidity  and  thus  with  the  refractive  index.  These  'artifi¬ 
cial'  features  can  be  turned  to  advantage  in  field  measurements  of  natural 
aerosol  since  they  provide  a  measure  of  the  refractive  index. 
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1  INTRODUCTION 

Pinnick  and  Auvermann'  have  discussed  the  response  characteristics  of  various 
Knollenberg  aerosol  counters  -  in  particular  the  active  scattering  aerosol  spectrometer 
(ASAS)  and  the  classical  scattering  aerosol  spectrometer  (CSAS) .  They  stress  that  the 
response  is  sensitive  to  the  aerosol  refractive  index.  Moreover,  even  with  a  given 
refractive  index,  spherical  particles  of  different  sizes  can  produce  the  same  response  in 
certain  regions;  these  are  described  as  regions  of  multi-valued  response,  although  it  is 
not  the  response  which  is  multi-valued  but  the  size  of  particle  to  be  inferred  from  a 
given  response. 

2  PRINCIPLES 

Both  counters  act  by  collecting  the  light  scattered  as  each  individual  particle 
passes  through  the  beam  from  a  He-Ne  laser.  The  CSAS  probe  is  conceptually  simpler  since 
the  scattering  volume  lies  outside  the  optical  cavity  of  the  laser.  As  shown  in  Fig  Jb, 
the  acceptance  solid  angle  of  the  collecting  optics  is  a  hollow  cone  between  half-angles 
of  4°  and  22°  off  the  forward  direction  of  the  laser  beam  -  the  missing  centre  of  the 
cone  being  determined  by  the  presence  of  a  dump  spot  to  stop  the  direct  beam  from  the 
laser.  In  the  ASAS  'active'  probe  (Fig  la)  the  dump  spot  is  replaced  by  a  plane  mirror 
so  that  the  scattering  volume  now  lies  within  the  optical  cavity  of  the  laser  (which  now 
has  a  Brewster  window).  The  geometry  of  the  collecting  optics  is  otherwise  identical  to 
the  CSAS,  although  the  particles  are  now  illuminated  by  light  travelling  in  both  direc¬ 
tions  with  consequent  interference  effects. 

Since  the  acceptance  solid  angle  of  the  collecting  optics  is  much  less  than  the 
tull  4tt,  it  is  clear  that  the  total  scattering  cross-section  is  not  strictly  applicable. 
Nevertheless  both  counters  collect  a  large  part  of  the  total  scattered  light  except  for 
particles  considerably  smaller  than  the  wavelength,  when  the  angular  distribution  tends 
to  the  ( 1  +  cos“  9)  form,  or  very  much  larger  than  the  wavelength  when  the  scattering 
becomes  concentrated  about  the  forward  direction.  We  shall  therefore  explain  the 
principles  underlying  the  effects  of  interest  by  considering  the  total  scattering.  The 
theory  is  conventionally  couched  in  terms  of  the  'scattering  efficiency  factor'  Qs(x,  m)  > 
where  x  =  2irr /\  ,  and  m  =  n^  -  in,  ,  where  n^  is  the  refractive  index  and  n2  is 

the  absorption.  By  definition  Q  (x,  m)  is  the  ratio  of  the  total  scattering  cross- 

.  .  S  .  2  . 
section  to  the  geometrical  cross-section  ^r  ,  te  the  total  scattering  cross-section 

2 

is  rrr  Q  (x,  m)  . 
s 

2 

Born  and  Wolf  state  that,  provided  f n j  -  1)  is  not  too  large  and  n,  is 
negligible,  Q  (x,  m)  is  approximately  a  universal  function  of  (n.  -  I ) x  with  its 

S  ] 

first  and  largest  peak  at  (n^  -  l)x  =  2  .  Fig  2  shows  the  total  scattering  efficiency 
factors  for  the  refractive  indices  Oj  =  1.3,  1.4,  1.5  and  1.6  ,  and  with  negligible 
absorption,  plotted  as  functions  of  the  dimensionless  argument  ( n f  —  1 ) x  .  The  func¬ 
tions  are  calculated  using  a  modified  form  of  the  standard  Mie  scatter  programmes  due 
3 

to  Dave  ;  since  the  programmes  require  that  n,  should  not  be  exactly  zero,  the  absorp- 
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tion  was  taken  as  n2  ”  10  .  Clearly  Fig  2  confirms  the  assertion^  that  the  gross 
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forms  of  the  functions  are  very  similar,  although  there  are  quite  marked  small-scale 

oscillations  in  the  individual  curves.  Fig  3  similarly  shows  the  total  scattering  cross- 
2 

section,  irr  Qg(x,  m)  ,  for  the  same  values  of  refractive  index,  now  plotted  against 

particle  radius  assuming  X  =  0.6328  ym.  The  horizontal  scale  is  chosen  to  include  ASAS 

range  0,  namely  0.3  to  1.5  pm  radius,  and  the  Y-axis  gives  the  cross-section  in  units  of 
2 

(pm)  on  a  log  scale.  All  the  cross-section  curves  show  a  plateau  or  even  a  slight  dip 
(confused  somewhat  by  the  small-scale  oscillations),  which  correspond  to  the  region  from 
the  first  gross  maximum  to  the  first  gross  minimum  of  Qg(x,  m)  in  Fig  2.  It  is  clear, 
however,  that  if  the  counter  were  detecting  the  total  scattered  energy,  then  a  range  of 
particle  size  in  the  region  of  the  plateau  or  dip  would  give  much  the  same  response. 
Moreover  the  particle  size  at  which  the  plateau  or  dip  occurs,  and  also  the  response, 
decrease  as  the  refractive  index  increases. 

It  must  be  emphasised  that  Figs  2  and  3  are  given  to  illustrate  the  principles  and 
are  not  strictly  applicable  to  the  Knollenberg  counters.  The  proper  calculation  of  the 
Knollenberg  counter  responses  is  a  much  longer  process  which  requires  the  calculation  of 
the  scattering  angular  distribution  followed  by  an  integration  between  the  limits  imposed 
by  the  geometry  of  the  detecting  optics.  The  results  of  such  exact  calculations  are 
given  by  Pinnick  and  Auvermann1  in  their  Figs  13  and  12  for  the  ASAS  and  CSAS  counters, 
respectively,  for  various  values  of  m  :  m  =  1.33  -  Oi  (water),  m  =  1.50  -  Oi  (ammonium 
sulphate),  m  *  1.50  -  0.00 5i  and  1.50  -  0.05i  (atmospheric  dust)  and  m  =  1.95  -  0.66i 
(carbon).  They  also  give  results  for  various  calibration  substances  (polystyrene,  glass, 
nigrosin  dye),  sodium  chloride  and  potassium  chlorate  in  their  Figs  5,  8  and  9  for  the 
ASAS  counter,  and  similarly  in  Figs  4,  6  and  7  for  the  CSAS  counter. 

The  exact  calculations  show  a  much  more  pronounced  dip  in  the  response  curves  for 
the  4°  to  22°  acceptance  cone  of  the  counters  than  is  found  for  the  total  scattered 
energy.  The  effect,  however,  is  the  same.  In  the  region  of  the  dip,  the  counter  necess¬ 
arily  lumps  together  particles  of  quite  distinct  sizes  because  they  give  the  same 
response.  This  'multi-valued'  region  is  followed  by  a  'single-valued'  region  where  the 
response  does  increase  monotonically  with  size  apart  from  the  small-scale  oscillations. 

The  effect  is  that  the  particle  size  distribution  inferred  from  measurements  contains  a 
'knee'  followed  by  a  sharp  drop.  As  has  been  emphasised  in  the  foregoing  discussion, 
and  also  in  Ref  1,  the  position  of  the  knee  should  depend  on  the  aerosol  refractive 
index. 

3  RESULTS 

We  have  made  various  sets  of  field  measurements  using  these  two  Knollenberg  aerosol 
counters  which  confirm  the  behaviour  described.  Here  we  use  data  taken  at  RAE  Farnborough 
from  25  September  to  3  October  1980,  as  the  UK  contribution  to  the  collaborative  aerosol 
measurement  progratmne  on  and  around  the  North  Sea  in  Autumn  1980  under  the  auspices  of 
NATO  AC/243  Panel  IV  R.S.G.8.  Our  active  probe,  the  ASAS- 300,  has  four  ranges  covering 
0.075  to  1.5  pm  nominal  radius,  while  our  classical  high-volume  probe,  the  CSAS-100-HV, 
has  only  two  ranges  covering  0.25  to  16  um  nominal  radius.  The  six  ranges  overlap 
considerably;  the  limits  will  be  sufficiently  clear  from  Fig  4  where  the  ranges,  each 


divided  into  15  nominally  equal  size  channels,  are  distinguished  by  different  symbols. 

To  ensure  proper  sampling,  the  probes  were  mounted  on  a  frame  which  was  controlled  to 

.  .  .  .  4 

point  into  the  wind  using  a  modified  commercial  aerial  rotator  system  slaved  to  a  wind 

vane.  The  whole  assembly  was  supported  on  a  tripod  with  the  probes  about  2  m  above  ground 

level  in  a  fairly  flat  grassy  area. 

Fig  4  shows  a  plot  on  a  log-log  scale  of  the  particle  size  distribution  n  =  dN/dr 
against  r  obtained  from  a  period  of  autorange  data,  in  which  the  probes  cycle  auto¬ 
matically  through  their  ranges,  on  29  September  1980.  This  example  was  deliberately 
chosen  as  the  autorange  period  which  shows  the  most  marked  knee  -  the  drop  occurs  almost 
entirely  between  channels  9  and  10  of  ASAS  range  0,  just  at  1  pm  nominal  size,  and 
exceeds  one  order  of  magnitude.  Where  they  overlap,  the  results  from  different  ranges, 
particularly  within  the  ASAS  probe,  are  in  excellent  agreement.  Also  CSAS  range  1 
(symbol  +)  agrees  well  with  ASAS  range  0  (symbol  *),  which  it  completely  encompasses, 
and  equally  shows  a  large  drop  in  imputed  number  density  near  1  pm  although  at  much  lower 
resolution.  Usually  the  knee  is  less  marked  than  in  Fig  4  because  the  drop  is  spread 
over  three  or  four  channels  perhaps  due  to  a  range  of  refractive  indices.  Even  for  a 
single  refractive  index,  however,  the  drop  should  be  spread  over  three  channels  if  the 
knee  is  really  in  the  middle  of  a  channel  rather  than  on  the  boundary  between  two 
channels . 

Usually  the  probes  were  both  operated  on  range  0,  giving  nearly  continuous  sampling 
of  ASAS  range  0  over  9  days.  To  get  adequate  statistics  the  data  was  then  averaged  hour 
by  hour.  As  an  automatic  computational  guide  to  the  position  of  the  knee,  we  simply 
found  the  ratios  of  the  counts  in  successive  channels  and  selected  the  largest  ratio. 

The  solid  line  at  the  top  of  Fig  5  shows  the  channel  number  preceding  the  knee  as  deter¬ 
mined  by  this  ratio  method  but  with  12  changes  made  after  visual  inspection  of  the  plots, 
the  original  values  being  shown  by  the  symbol  in  six  instances  the  largest  ratio 

is  quite  unconnected  with  the  knee.  Three  times  channel  14  is  chosen  by  statistical 
vagary,  twice  channel  1  is  chosen  because  there  is  a  steep  gradient  at  the  bottom  end 
of  range  0  and  the  drop  is  spread  over  a  few  channels,  and  once  channel  4  is  chosen  by 
a  confcination  of  these  effects.  There  are  also  six  minor  changes  in  the  last  three  days, 
where  the  drop  is  spread  over  a  few  channels  and  the  largest  drop  is  not  the  first  one. 

The  ratio  method  is  thus  not  entirely  successful,  but  no  doubt  a  better  algorithm  could 
be  devised.  It  is  clear,  however,  that  the  position  of  the  knee  varies  systematically 
with  time. 

Since  the  position  of  the  knee  ought  to  depend  on  the  refractive  index,  it  should 
therefore  be  correlated  with  relative  humidity.  Air  temperature  and  dew  point,  and  thus 
RH,  were  measured  using  an  EG4G  Model  220  Dew  Point  Hygrometer,  and  automatically 
recorded  every  100  s.  The  lower  curves  in  Fig  5  show  the  average  RH  over  each  hour 
together  with  the  extreme  values  during  the  hour.  Clearly  the  position  of  the  knee  does 
indeed  correlate  well  with  RH.  Table  1  gives  the  correlation  coefficients  (i)  using  the 
channel  number  given  by  the  ratio  method,  (ii)  after  making  the  obvious  changes,  and 
(iii)  after  making  all  changes.  Values  are  given  for  the  data  overall,  and  also  split 
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up  into  groups  of  3  days.  The  obvious  changes  from  channel  1  (twice)  and  channel  4  (once) 
have  the  greatest  effect  on  the  correlation,  while  the  minor  changes  have  almost  no  effect 
Both  the  second  and  third  groups  of  three  days  contain  periods  when  both  the  external 
laser  mirror  in  the  ASAS  probe  and  the  misting  mirror  in  the  GG&G  dew  point  sensor 
obviously  warranted  cleaning  and  duly  were  cleaned  (see  gaps  in  Fig  5);  this  must  partly 
explain  the  lower  correlations  in  Table  1.  The  correlation  is  extremely  high,  however, 
given  that  the  position  of  the  knee  is  quantised,  that  Oj  is  unlikely  to  depend 
linearly  on  RH,  and  that  the  presence  of  different  substances  could  give  distinctly 
different  values  of  n|  • 

Table  1 


Correlation  coefficient  between  the  position  (channel  nunfcer)  of  the  'knee' 
on  ASAS  range  0  and  the  relative  humidity  based  on  measurements 
at  RAE,  Farnborough  in  September/October  1980.  The  number  of 
changes  is  given  in  brackets 


25-27  September 

28-30  Septenber 

1-3  October 

Overall 

Given  by  ratio  method 

0.627 

0.454 

0.807 

0.610 

After  obvious  changes 

0.930  (2) 

0.755  (2) 

0.823  (2) 

0.821  (6) 

After  all  changes 

0.930  (2) 

1 

0.755  (2) 

0.852  (8) 

0.844  (12) 

Nunber  of  points 

57 

63 

58 

178 

As  shown  in  Fig  5,  the  lowest  position  for  the  knee  is  after  channel  7  on 

25  September  when  RHav  is  55.3%  and  on  26  September  when  RHav  is  61.6%.  It  clearly 

lies  after  channel  14  twice  on  2  October  when  RH  is  99.2%  and  99.3%.  In  their 

1  .  av 

Fig  13,  Pinmck  and  Auvermann  give  the  results  of  exact  calculations  for  the  response 
of  the  ASAS  counter  to  particles  with  m  =  1.33  -  Oi  (water)  and  m  =  1.50  -  Oi  (specifi¬ 
cally  ammonium  sulphate  but  reasonably  typical  of  other  hygroscopic  substances)  assuming 
that  the  pulse  height  discriminator  levels  are  set  as  recommended  by  the  manufacturer . 
From  their  Fig  13,  the  knee  should  appear  after  channel  15  (for  m  =  1.33  -  Oi)  and 
channel  7  (for  m  *  1.50  -  Oi)  respectively  on  ASAS  range  0.  If  our  discriminator 
settings  were  also  as  recommended  by  the  manufacturer,  our  highest  and  lowest  knees 
correspond  to  refractive  indices  of  about  1.35  and  1.50  respectively,  by  interpolation 
from  the  curves  in  Ref  I. 

In  Fig  4  further  oscillations  are  evident  on  CSAS  range  1  with  peaks  near  2  urn  and 
3  pm  radius  which  correspond  to  the  further  oscillations  in  the  scattering  efficiency 
factor  (see  the  topmost  curve  in  Fig  6)  which  the  counter  cannot  resolve  on  range  0. 
Specifically  the  peaks  in  Fig  4  should  coincide  with  the  portion  between  each  maximum 
in  Fig  6  and  the  succeeding  minimum.  These  oscillations  could  likewise  be  used  to  find 
the  refractive  index  for  the  relevant  sizes  of  particle,  but  no  analysis  is  attempted 
here  since  the  range  is  sampled  only  occasionally  in  the  present  data.  The  oscillations 
on  CSAS  range  1  could  also  give  some  indication  of  the  absorption  since  the  corresponding 


041 


peaks  in  the  scattering  efficiency  factor  decrease  with  increasing  absorption  and  have 
effectively  vanished  when  =  0.1  as  shown  in  Fig  6.  We  have  verified  that  the 
oscillations  are  not  present  if  the  counter  is  operated  indoors  in  a  low  RH  atmosphere 
when  the  aerosol  particles  will  be  solid,  dusty  and  absorbing. 

4  CONCLUSIONS 

We  therefore  conclude  that  the  distortions  imposed  on  the  inferred  particle  size 
distribution  by  the  characteristics  of  the  counter  responses  are  really  an  advantage  in 
field  measurements.  The  position  of  the  knee  on  ASAS  range  0  provides  a  good  measurement 
of  the  refractive  index  (at  the  He-Ne  wavelength)  which  could  be  put  on  an  even  sounder 
footing  by  regular  electronic  calibration  or  monitoring  of  the  discriminator  settings, 
etc.  The  distortions  also  provide  a  sensitive  test  that  the  probes  are  operating 
properly,  and  almost  certainly  could  be  adequately  corrected  for  in  the  data  analysis. 

In  any  case  the  distortions  arise  from  the  form  of  the  scattering  efficiency  factor  for 
particular  ranges  of  particle  size,  and  in  practical  terms  seem  unavoidable  in  an  optical 
device. 

Pinnick  and  Auvermann'  have  suggested  that  channels  should  be  aggregated  together 
in  preset  groupings  -  no  doubt  a  reasonable  procedure  if  the  refractive  index  is  known 
and  remains  constant.  With  field  measurements,  however,  any  grouping  of  the  channels 
should  be  done  only  after  the  position  of  the  knee  has  been  determined  either  from  the 
size  distribution  or  even  by  inferring  the  refractive  index  from  the  relative  humidity. 

If  the  results  from  the  aerosol  counters  can  be  shown  to  be  reliable,  and  if  the 
sampling  volume  and  discriminator  settings  are  known,  the  results  could  be  used  to  find 
the  extinctions.  Knowing  the  refractive  index  the  real  sizes  of  the  particles  could  be 
found.  At  the  same  time  the  ’knee’,  and  also  the  further  oscillations  if  necessary, 
could  be  removed  by  curve-fitting.  The  extinctions  at  various  different  wavelengths 
could  then  be  calculated.  This  Report  is  a  step  in  that  process. 
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Fig  4  Particle  size  distribution  from  a  period  of  autorange  data  on  29  September  1980, 
at  RAE,  Farnborough.  Abscissa  is  the  radius  of  the  channel  centre,  and  the 
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ranges  0  and  1  of  the  CSAS  probe  respectively 
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Fig  5  Time  plot  of  the  position  of  the  knee  on  ASAS  range  0  and  the  relative 
humidity  (min,  av,  max)  during  hour.  All  data  averaged  over  1  hour 
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Fig  6  Total  scattering  efficiency  factors  with  increasing  absorption  (*-  0, 0.03, 0.1) 
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